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A nonequilibrium master equation analysis for TV interacting bosons, with Bogoliubov quasiparti- 
cles as the reservoir is presented. The analysis is based on a simplified Hamiltonian. The steady state 
solution yields the equilibrium density matrix. The results are in good agreement with and extend 
our previous rigorous canonical ensemble equilibrium statistical treatment leading to a quantum 
theory of the atom laser. 



I. INTRODUCTION 



We here analyze the condensate fluctuations of TV in- 
teracting bosons via a master equation analysis based on 
a simple truncated Hamiltonian. The strategy is to re- 
gard excited states of the Bose g 0jS, 0jS well as collective 
excitations, as reservoir variables which we ultimately 
trace over, in order to obtain a master equation for the 
condensate. We find that the steady state solution of the 
resulting kinetic equations provides a good description 
of Bose-Einstein condensate (BEC) Fig. ^ shows the 
condensate fraction as a function of temperature and is 
one of our main results. This is in accord with the wit 
and wisdom of Wigner: 

With classical thermodynamics, one can cal- 
culate almost everything crudely; with kinetic 
theory, one can calculate fewer things, but 
more accurately; and with statistical mechan- 
ics, one can calculate almost nothing exactly. 
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In particular, we discuss several subtle issues concern- 
ing the statistics of condensate atoms in the region of 
T ~ T c . This is analogous to studying the photon statis- 
tics of the laser in the passage from below to above 
threshold. In fact, BEC is often referred to as an atom 
laser, and the question of the atom statistics of the meso- 
scopic BEC (containing a few 10's or 100's of atoms), is 
the subject of this paper. The answer is given by the di- 
agonal elements of the BEC density matrix, p„ 0i „ , where 
no is the number of atoms in the BEC ground state. In 
the present case of a gas of N weakly interacting bosons, 
in one limit we find 



FIG. 1: Solid lines are the present work results showing the 
mean number of condensate particles as a function of temper- 
ature for N = 200 particles in a box calculated via the solution 
of the condensate master equation for an interacting (Bogoli- 
ubov) gas with an 1,/3 =0.1 and an ideal (Bose) gas are 1 ' 3 = 0. 
Dots are the exact numerical results obtained in the canonical 
ensemble for the ideal Bose gas Q. Squares are the results 
of CNB3 (please note that repulsive interaction increases T c ). 
Dashed line is a plot of 7V[1 - (T/T c ) 3/2 ] which is valid for 
the ideal gas in the thermodynamic limit. The temperature 
is normalized by the ideal gas thermodynamic critical tem- 
perature in the box T c = 2irh 2 n 2/3 /k B M((3/2) 2/s , where M 
is the particle mass. 
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where Zn is a normalization (partition function) factor 
expressed in terms of the incomplete gamma function 
T(N + 1, Hi); and Hi is a simple heating coefficient gov- 
erning the rate of removal of atoms from the ground state 
due to scattering by Bogoliubov quasiparticles, which is 



given by 



Hi = 



(2) 
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In the above the usual Bogoliubov amplitudes are 

^ k = ^l £k " W"— J' (3) 

where M is the atomic mass, V is the condensate volume 
and [/k is the atom-atom scattering energy; the quasi- 
particle energy e k is given by 

l(W n U k \ 2 fn U k \ 2 

The simple expression, Eq. ft}, is valid in the limit 
in which, for the case of an ideal gas in a box, Tij = 
N (T/T c ) 3 ^ 2 where T is the temperature of the gas and 
T c is the usual critical temperature. This was the main 
result of the first paper on the Condensation of N Bosons 
(CNB1) in this series Q, in which we focused on the 
statistics of N noninteracting (ideal) bosons. In the sec- 
ond paper 4] (CNB2) the approximations were improved, 
as discussed in section III below. In CNB1,2 we solved for 
the condensate statistics of an ideal Bose gas and found, 
surprisingly, that the fluctuations are not Gaussian even 
in the thermodynamic limit. The treatment of CNB1,2 
is patterned after the quantum theory of the laser, i.e., 
is based on a master equation for the condensate density 
matrix. As such, it has the useful (and intriguing) feature 
of working equally well at all temperatures from T = 
to T c and above. Atom-atom interactions are, however, 
ignored in those papers. 

Paper CNB3 .5] includes interactions. The analysis 
of CNB3 is based on the Girardeau- Arnowitt adaptation 
of Bogoliubov's approach to BEC. In CNB3, we give a 
"rigorous" statistical mechanical treatment of the fluc- 
tuations, valid provided the average number of conden- 
sate particles is much larger than its variance; that is 
valid for temperatures not too near T c . This condition 
implies that the excited atom numbers n& fluctuate in- 
dependently, < rifcn m >= Ufcn m , k ^ m, by exchanging 
particles with the condensate reservoir. For N = 200 
CNB3 is valid from T = to T w 0.7T C . However, near 
and above T c the correlations become substantial which 
yields failure of such a treatment. In a broad tempera- 
ture range, interactions increase the condensate number 
(no), and the fluctuations are found to be one half those 
of the ideal Bose gas. This remarkable fact is due to the 

(jt — kj pairing of atoms. 

In the present paper, we utilize a master equation anal- 
ysis for the interacting Bose gas based on the simple trun- 
cated Hamiltonian given by Eq. 114fl . In particular we 
derive the average number of particles in the lowest en- 
ergy level, (no) and the fluctuation about the average 
An 2 = ((no - (no)) 2 . These are plotted in Figs. Q]and[2] 
and constitute our main results. 

We note that the present analysis yields an expression 
for (no), in good agreement with the rigorous results of 
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FIG. 2: Variance Auq of the condensate particle number as a 
function of temperature for N = 200 particles in a box calcu- 
lated via the solution of the condensate master equation for 
an 1 / 3 = and 0.1 (solid lines). Dots are obtained by exact 
numerical calculations for the ideal Bose gas in the canonical 
ensemble Q|. Squares are the result of CNB3. For an 1 /' 3 = 
0.1 the solid line is obtained from the master equation by a 
shift of Hi to match the CNB3 value at T = (see Table 1 
column three). 

CNB3, in the temperature region where that analysis is 
valid; and extends those results to higher temperature up 
to T c and beyond, as per Fig. ^ 

We emphasize that the problem of the interacting Bose 
gas near T c is a difficult one. It is therefore worthwhile 
to investigate simple models designed to capture the es- 
sential physics. To that end, the next section is based on 
a simple "toy" model which works surprisingly well. En- 
couraged by the results of that analysis we then present 
a physically motivated model yielding a similar Hamil- 
tonian. The analysis based on a more general treatment 
will be published elsewhere. 

In the next section we derive the relevant master equa- 
tion which we solve in section III to obtain the condensate 
statistics. In section IV we discuss the results and make 
contact with previous work. 



II. MASTER EQUATION FOR WEAKLY 
INTERACTING BOSE GAS 

Having set the stage, we proceed to sketch a simple 
approach extending our master equation analysis to the 
case of the interacting Bose gas. The present treatment 
is based on simple extension of CNB1,2 via a "toy" model 
and then a more realistic (but still heuristic) formulation. 
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However the results are in good agreement with "rigor- 
ous" results of CNB3 for temperatures where that theory 
is valid and provide answers in the vicinity of T c where 
it does breakdown. 

In CNB1,2 the basic model interaction Hamiltonian 
describing the statistical dynamics of the condensate was 
taken to be 



V = 



5k a bl a k 



adj. 



(5) 



where o,k and 6k are atom and phonon annihilation op- 
erators, likewise do is the condensate annihilation oper- 
ator and is the corresponding coupling strength for 
the collision of a ground state atom and an atom having 
momentum k scattering into the BEC. This models the 
dilute gas BEC experiments of Reppy and coworkers @ . 
Moreover the results for the condensate particle number, 
obtained by the master equation approach of CNB2 and 
exact numerical simulation for the ideal gas in the canon- 
ical ensemble |2j , are in excellent agreement as indicated 
in Fig. [2 



A. Toy Model for an Interacting Bose Gas 

As a first step toward obtaining no and Ano for an 
interacting Bose gas, we derive a master equation for 
Pn ,n by simply replacing the phonon operators 6k and 
6 t in Eq. (JSJ by operators 8k and B k related to 6k and 
6 k by the usual Bogoliubov transformation 



fi k = Uk6k + «k&Li 
fi k = u k 6 w + Wk6-i 



(6) 



where Uk and «k are given in Eq. J2J, < 6 k 6k >= 
l/(e' 3ek — 1), e k is the Bogoliubov quasiparticle energy 
given by Eq. @. 

The present master equation for the condensate is now 
obtained via the phenomcnological interaction Hamilto- 
nian 



V = Y #k «k fik «o + ad J 



(7) 



which adds (removes) atoms from the condensate with 
the annihilation (creation) of an excited atom and the 
emission (absorption) of a quasiparticle. 

We recall that the simple ideal gas master equation for 
the probability P na of finding no atoms in the N atom 
gas ground state was found to be 3] 



where k is an uninteresting overall rate factor and the 
phonon heating coefficient was given by 
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(9) 



where (3 = l/k B T. 

In the present toy model all we do is replace the phonon 
operators 6k and 6 k in Eq. (0 by the quasiparticle op- 
erators fik and fi k of Eq. ©, and find the new heating 
coefficient 
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E 
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(10) 



The results of this simple analysis are summarized in 
Figs. and Encouraged by these results we next 
present a simple argument which puts the phenomeno- 
logical analysis of this subsection on a more physically 
motivated footing. 



B. A Physical Argument Based on 
Bogoliubov-Girardeau-Arnowitt Formalism 

We work with the particle number conserving creation 
and annihilation operators of Girardeau and Arnowitt Q 

/3k = 4A), /3k = /?o«k, $ = (l+n )- 1/2 a . (11) 

These operators describe canonical-ensemble quasiparti- 
cles which obey the Bose canonical commutation rela- 
tions 



0k, ft 



?k,k' 



(12) 



see CNB3. The canonical-ensemble quasiparticles /3k = 
/?Qcik describe transitions between ground (k = 0) and 
excited (k ^ 0) states. 

Consider next the conventional atom-atom interaction 
Hamiltonian as sketched in Fig. [3] Thus we may write 

Vi = ^ Uv\a\_^alaoa\ = ^ f7kidj_ k /3k"i (13) 



k.l 



k.l 



where we have introduced the canonical-ensemble quasi- 
particles a k eto ~ V^oS^l + "-o) _1 ^ 2 ao = V™oftt.> and 
the notation — ^/noUid with Uk\ being the scattering 
matrix element. Finally we take 1 = k (leading term) 
and thus write the interaction Hamiltonian 



V r = ^C/ap k a k + adj. 



(14) 



-P na = -(N - n ){n Q + l)P no + (N-n + l)n P (lo _i- 

K 



Hn P no +H(n + l)P, 



n + li 



(8) 



It is important to note that this Hamiltonian involves 
two different kinds of operators. The operators a Q and ak 
are single particle operators which create and annihilate 
particles. The /3 k operator, by construction conserves 
particle number as befits a phonon like quasiparticle. 
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average ^/3 k /3 k y is replaced by an average number of ther- 
mal phonons in our earlier work. The physics behind i|15|) 
is best understood by going to a "low temperature limit" 

such that (plfik ) <C 1 and we may write Eq. i|16[l as 



(18a) 



'fi k K ^o = E ( n ^n„ =N-n , 

and (^k)„ C 1 so that we may write Eq. I|17|) as 



FIG. 3: Two particle scattering process involving the con- 
densate, and the same process in terms of canonical ensemble 
quasiparticles. 



C. The Master Equation 

The present master equation for the condensate is now 
derived by focusing on the truncated scattering Hamil- 
tonian Eq. (|14l) . which adds (removes) atoms from the 
condensate with the annihilation (creation) of an excited 
atom and the emission (absorption) of a quasiparticle. 
As in CNB1 and CNB2 we treat the phonon like quasi- 
particle states as reservoir states for the condensate. 

The dynamical evolution of the condensate thus in- 
volves three components: the ground state (conden- 
sate) atoms, excited (noncondensate) and the phonon 
like canonical ensemble excitations, /3k and /3 k . We pro- 
ceed to treat the /3 k quasiparticles as a thermal reservoir 
which we trace over to obtain the density matrix equa- 
tion for the ground state, as in CNB1,2. Some details 
of the master equation analysis are given in Appendix 
A. As will be discussed further elsewhere, the equation 
of motion for the probability to find no particles in the 
condensate p no ,n = Pn is now given by 



1 



P no = -K no (n + 1)P„ + K no ^in P no ^ 



H no n a P na + H no+1 (n Q + l)P na+1 . 



(15) 



As in CNB1,2, we have divided the physical process into 
two kinds of terms: the K no and K no -i terms describ- 
ing the cooling of the gas which increases the condensate 
number and heating terms H no and H no+ i which de- 
crease it. The constant k is an uninteresting overall rate 
factor. The cooling and heating coefficients are given by 



K no = («*>no f 1 + (®A 



k^O 



H no = E(^k)(l + (n k ) no ) 

k^O 



(16) 
(17) 



Equation i|15|) was obtained previously in a similar form 
for the ideal gas; the essential difference being that the 



Hr* = Hi. 

k^O 



(18b) 



We proceed by writing the /3k and /3 k operators in 
terms of the Bogoliubov quasiparticle operators 



Ac = tt k b k + w k&_ k 
0- k = "" k b^ k + v k b k 

where, as usual, [b k , bt,] — <5 k , k ', hence 



(19) 
(20) 



k \ u k u k/ T L' k ^u_ k u_ k 
«k-^k 



e 0e k _ I 



(21) 



In the preceding we have taken the quasi-particle to be 
in thermal equilibrium such that 



b k b k 
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s /3c k — 1 ' 



(22) 



Finally, we note that in the useful notation of CNB3 the 
Hi coefficient found from Eqs. (|18 h>|) and (|2"T)> reads 



1 



1 



Hl 2^ o Vz(A k )-l ' z(-A±)-lJ> 



where 



z(A h 



A^I T - 1 



(23) 



(24) 



and Ak is given by Eq. J3J). 

In the low temperature limit the cooling and heating 
coefficients are given by Eqs. Ijl8a|l and l|18b(l . However, 
as was shown in CNB2, the cross coefficient 



E ^)n (fiA 



k^O 



(25) 



as it appears in Eqs. (|16fl and (|17|l is necessary in order 
to obtain an accurate description of the condensate. As 
was shown in CNB2 and as sketched in the appendix B, 
for the ideal gas we have 

E (^)no (^k) = (N- n ) V (26) 

k^O 
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where the phonon average number (fe^k 



1) 1 = 77k and the CNB2 cross coupling coefficient n 
is defined by 



where 



k^O 



(27) 



(28) 



For the interacting gas the heating coefficient Hi and 
the cross coupling coefficient can now be written as 



k^O 



e /3e k _ 1 



2 1 — 



Z>k (29) 

k^O 



and 

E ^ 

k^O 

where 



ftA) = E (^)n„^ = (N-n )m, (30) 

k^O 



k^O 



(31) 



Written in terms of 7i/ and rjj the master equation for 
the interacting Bose gas is finally 

-P„ = -(JV - no)(l + fj/)(no + 1)P« + 



[iV-(no-l)](l + r?/)nQP flo _i-[Wi + (Af-n )??j]noP„o + 



[Wj + (JV - (n + l))n/](n + l)P„ +i 



(32) 



III. CONDENSATE STATISTICS VIA STEADY 
STATE SOLUTION TO MASTER EQUATION 

The BEC atom statistical distribution is now obtained 
by solving the master equation l|32|) . The resulting dis- 
tribution is given by essentially the same expression as 
was obtained for the ideal Bose gas in CNB2, namely 



p — 

r n — 



1 (N-no + Hi/vi-lV- ( m 



Z N (Hi/rn-l)l(N-no)\ \l + m 



N-n Q 



(33) 



JV 

^=E 



N - n + Hj/m - 1\ / m 
N-n J \l + m 



N-tiq 



where Hi and r\i for the interacting gas are given in Table 
1. 



In the low temperature limit we set r\i to zero and 
obtain the distribution given by Eq. (J*[J. The average 
condensate particle number from Eq. Q is 



(n ) - N - Hi + H^ +1 /Z N N\ 



(35) 



This is in good agreement with the average value ob- 
tained in CNB3 namely (no) = N — Hi valid for T not 
in the vicinity of T c . The extra term Hf +1 / 'ZjyNl as 
given in Eq. (|35|l removes the unphysical cusp at T = T c 
and extends the treatment of CNB3 through the critical 
temperature. Equation Q and the associated (no) is es- 
sentially the extension of the results of CNB1 to include 
atom-atom interaction. 

In general the r\i parameter is important and the qua- 
sithermal (no) obtained from Eq. (|33|l is 

(n ) =N-Hi + P {mN + Hi) 1 



it is plotted in Fig. ^ 
The squared variance 

Ann = 



n 



(36) 



can be calculated analytically from Eq. I|33|) to find 
An 2 = (1 + m )Hi - P (viN + Hi){N -Hi + 1 + m) 



where 



^0 



P*{ m N + Hif 



1 (N + Hi/ m -iy. 



Z N N\(Hi/r]i-l)\ \l + Vi 



Vi 



N 



(37) 



(38) 



is the probability that there are no atoms in the conden- 
sate. 

If the temperature is not too close to the critical tem- 
perature, only the first term in Eq. (|37(l remains, result- 
ing in 



An 2 Q ~ (1 + m )Hi 



k 



(mi + <n k ) 



(39) 



We plot the variance in Fig. |2Jfor the case of 200 atoms 
in a box and the values of gas parameter an 1 / 3 = 0, 0.1. 
Here a = MU / Airti 2 is the s-wave scattering length, M is 
the atom mass, U ~ [/ki is the scattering matrix element 
in Eq. i|13|l and n is the particle density in the box. 
Squares show the result obtained from CNB3, that is 



((n - n ) 2 ) = - 



2 V 

k#0 



1 



1 



(z(Ak)-l) 2 (z(-Ak)-l) 5 



z(A k ) - 1 z(-A k ) - 1 
which is in agreement with Ref. |8j]. 



(40) 
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low temperature 
approximation 


quasithermal 
approximation 


quasithermal 
approximation 
improved by CNB3 


Hi 













EfcVo ^ 


E M0 ^ 



TABLE I: Expressions for Hi and 77/ in three limits of opera- 
tion corresponding to the three previous treatments extended 
to the present model. The average particle number in the 

excited level is given by 9^= 



Table 1 shows Hi and r\i in different approximations. 
The interacting Bose gas values of Hi and nj are given 
in column one in the CNB1 low temperature limit in 
which rji = 0. The second column gives the heating Hi 
and cross coupling r\i coefficients according to the Eqs. 
(|29J) and (|31|) . However, we have specific information 
from CNB3 such that we know Ano etc. exactly at, say, 
T = 0. We use that information to improve the present 
treatment of Ano by shifting Hi by a constant value A 
so that we match Ano of CNB3 at T = 0. 

Figure 21 compares the variance Ano of the condensate 
particle number calculated for an 1 / 3 = 0.1 in three dif- 
ferent approximations. The dash-dot line is obtained in 
the low temperature limit (Table 1 column one). The 
dashed line is the quasithermal approximation (Table 1 
column two) . The solid curve was obtained from Hi and 
rji of Table 1 column three. 
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FIG. 4: Variance Ano = \J (n^) — (no) 2 of the condensate 
particle number as a function of temperature for N = 200 
particles in a box calculated via the solution of the conden- 
sate master equation for an 1 ^ 3 = 0.1 in different approxima- 
tions. Dash-dot line is obtained in the low temperature limit. 
Dashed line is the quasithermal approximation. Squares show 
the result of Ref. y| and CNB3. Solid line is calculated in 
the quasithermal approximation by a shift of Hi to match the 
CNB3 value at T = 0. 



IV. DISCUSSION AND SUMMARY 
A. Mean number of particles in the condensate 

The calculation of (no) for the interacting "Bose- 
Bogoliubov" gas is in good agreement with and extends 
the results of CNB3. For example, the Uhlenbeck cusp 
dilemma is resolved for the interacting gas just as it 
was for the ideal gas. The master equation approach 
gives an excellent treatment of the temperature depen- 
dence average number of atoms in the condensate. The 
agreement with the rigorous results of CNB3 is grati- 
fying. Note in particular the ground state depletion at 
T = 0, due to atom-atom interactions, is handled very 
well by the present approach. Likewise the increase in 
the number of atoms in the ground state, for temper- 
ature 0.3 < T/T c < 1.0, is modelled very well by the 
master equation approach. For temperatures well above 
T c (T/T c > 1.4) the ideal gas treatment is valid and 
agrees with our result. Furthermore, and most signifi- 
cantly from the vantage point of this paper, the master 
equation approach works well for the temperature region 
near T c where the other treatments fail. 

As one can see from Fig. ^ repulsive interaction be- 
tween Bose particles stimulates BEC and yields an in- 



crease in no at intermediate temperature as compared to 
the ideal gas 01 • The reason is energetic: bosons in dif- 
ferent but mutually overlapping states interact stronger 
than when they are in the same state. For example, when 
two particles are in the same state ^(ri, r 2 ) = x(ri)x(r 2 ) 
and interact via a potential V = g5{r\ — r 2 ), the interac- 
tion energy is 

E- mt = I dr 1 /dr 2 **(r 1 ,r 2 )y*(r 1 ,r 2 ) = 



9 I dri|x(ri)| 4 . 

However when two particles are in different states such 
that *(ri,r 2 ) = [^(ri)x(r 2 )+^(r 2 )x(ri)]/\/2, the inter- 
action energy is 

E int = 2g J dnl^n)! 2 ^^)! 2 . 

Then since |<^(ri)| 2 ~ |x( r i)| 2 we see that two bosons in 
the same state is the lowest energy configuration. This 
favors multiple occupation of a single one-particle state. 
Such an effect is sometimes called an attraction in mo- 
mentum space |l l| . One can see from Fig. ^that repul- 
sive interaction also increases T c . 
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B. Fluctuations Ano 

As Einstein taught us long ago, fluctuation phenom- 
ena contains much more physics then mean values do. 
He used the difference between the fluctuation proper- 
ties of waves and particles to show the particle side of 
the photon an d the wave side of matter Note 
that the latter was well in advance of the Schrodinger 
equation, and provided support for the wave nature of 
matter. 



The extension of the ideal gas work of CNB2, as per 
column 2 of Table 1 to the present problem yields the 
dashed curve which is better but not good for T = 0. 
However, CNB3 provides a good description of the cen- 
tral moments ((no — uq)" 1 ) for low temperatures. Thus 
we can improve our description of the fluctuations by us- 
ing the known low temperature behavior from CNB3 to 
improve the heating coefficients as in column 3 of Table 1. 
The utilization of the low temperature results of CNB3 to 
improve on K no and H na yields an excellent description 
of the central moments, as will be shown elsewhere. 
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FIG. 5: Variance Ano = \J (n 2 ,) — (no) 2 of the condensate 
particle number as a function of temperature for an ideal Bose 
gas of N = 200 atoms in an isotropic harmonic trap. Solid line 
is a solution of the condensate master equation. Dots are the 
exact numerical results obtained in the canonical ensemble for 
the ideal Bose gas 0. Squares are result of CNB3. Dashed 



line is a plot of Ano = \J ^(3^ which is obtained in the 

thermodynamic limit 14]. Dash-dot line is a plot of Ano = 
N — (no) = N (jf^ , which is proposed by D. ter Haar [Tr| | 
in the low temperature regime (adapted to a harmonic trap). 
This had the correct zero limit as T — » 0, but is not right for 
higher temperatures. 

Likewise fluctuations are central to our investigation. 
We find that even in the ideal Bose gas the fluctuation 
physics is quite rich, see Fig. The atomic statistical 
distribution, and its first two central moments, as plotted 
in Fig. (for no) and Fig. [2] (for Ano) make our point. 
The low temperature approximation used to get Eq. l[TJl 
works quantitatively well for no. However, the treatment 
of fluctuations in this approximation is only qualitative 
(see Fig. 0J. To get a good description of fluctuations 
we must extend the improved approach of CNB2, with 
appropriate modifications, to the case of the interacting 
gas. 



Connection with previous studies and future 
work 



The experiments in dilute Bose gases of rubidium, 
lithium and sodium have reopened the old BEC ques- 
tions, e.g., the Uhlenbeck cusp dilemma Q and the ques- 
tion of condensate fluctuations. See for example the re- 
cent review |l6j and references therein. Likewise the con- 
densate time development [l7^ is an exciting frontier. 

In earlier work Q, we found that there is a use- 
ful connection between BEC of an ideal Bose gas [Isj . 
and the quantum theory of the laser 0, |2£j ■ We recall 
that the saturation nonlinearity in the radiation matter 
interaction is essential for laser coherence pi). 

Indeed, the coherence generating nonlinearity in the 
case of the ideal Bose gas is the particle number con- 
straint which provides the essential nonlinearity in BEC. 
However, in the case of the interacting gas, the problem 
of fluctuation is rather subtle and requires a more refined 
approach. Although useful papers have been published 
dealing with various limiting cases H, l^, so far there 
has been no treatment of this problem valid at all tem- 
peratures. Apart from general theoretical interest, con- 
densate fluctuations can be measured, in principle, by 
means of a scattering of series of short pulses [23, see 
also 0. 

The problem of N ideal bosons in a 3D harmonic po- 
tential coupled to a thermal reservoir turns out to be 
surprisingly rich. The N particle constraint is included 
naturally in the present formulation and introduces an 
essential nonlinearity yielding accurate results as shown 
in Fig. [T]and|21 

Given the utility of the quantum theory of the laser- 
master equation approach in the ideal Bose gas problem 
it is natural to look at the interacting gas from a quantum 
optical perspective, and that was the vantage of Ref. [j| . 
Other fascinating studies in this regard include the paper 
by Lewenstein and You on quantum phase diffusion [25j , 
Graham j2(| who applied quantum noise analysis to the 
BEC lincwidth question and Wiseman and Thomson on 
reducing the linewidth of an atom laser by feedback [27| • 
The time dependent master equation analysis of Gardiner 
and Zoller is presented in a heroic series of strong papers 
28]. The paper of Parkins and Walls [U on the BEC is 
another important work in this field. Reference to many 
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other useful and insightful papers is given in our recent 
review on fluctuations in the BEC [Ig . 

One basic difference between the classic papers men- 
tioned above is treatment of the BEC of a mesoscopic 
gas in the critical region near T c . The present master 
equation approach is, to our knowledge, unique in that 
it treats BEC analytically at all temperatures. 

In the next paper (CNB5) we shall follow the lead of 
Section ITTT1 see Table I and Fig. wherein we use the 
statistical mechanical results of CNB3 to find a good 
(semi-phenomenological) description of all moments at 
arbitrary temperatures. Then in CNB6 we shall present 
a detailed analysis including the off-diagonal elements 
(e.g. apajpo). Finally we will return to the study of 
various forms of the dynamics (Hamiltonian) in CNB7, 
and show that certain formulations do not involve off- 
diagonal elements and give a good account of the physics 
via a simple treatment. 
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APPENDIX A: THE MASTER EQUATION 



In CNB1,2 the equation of motion for the total 
"gas+reservoir" system density matrix in the interaction 
representation was given by 

Ptotal (t) = -i[V(t) ) Ptotal (t)]/h. (Al) 

Integrating for ptotai, inserting it back into the equation 
of motion for p(£)totai, and tracing over the reservoir, we 
obtain a useful equation of motion for the density matrix 
of the Bose-gas subsystem 

= -Jpf dt'Tr Ies [V(t), [V(t'),p tom (t')}], (A2) 

where Tr res denotes the trace over the /3k degrees of free- 
dom. 

We proceeded by assuming that the phonon reservoir 
(in the spirit of the experiments of Reppy and coworkers 
6]) remains unchanged during the interaction with the 
Bose gas. As discussed in detail in CNB1,2, the density 
operator for the total system "gas+reservoir" can then 
be factored, i.e., p to tai(£') ~ p{t') (E> Pros, where p ros is the 
equilibrium density matrix of the reservoir. If the spec- 
trum is smooth, we are justified in making the Markov 
approximation, viz. p(t') — > p(t). 

In the present analysis we do not invoke an external 
(phonon) reservoir. Instead we regard the excited states 



of our weakly interacting gas to be reservoir-like. In 
this regard the excited states are nearly analogous to the 
atoms in the quantum theory of the laser . There we 
calculate the change in the laser radiation density matrix 
due to one atom and multiply by the rate of such atomic 
interaction to obtain a coarse-grained equation of motion 
for the laser density matrix. 

Following the laser approach explicitly (in CNB1,2 it 
was only implicitly) we first calculate Sp(ao,a,Q,t) due 
to a collision between two atoms, which adds (or re- 
moves) atoms to the ground state and then multiply by 
the rate of such collisions to find our master equation for 
p(a,Q, <2q, t). The small change 5p is given by 

t + T c t + t' 

Sp(a ,alt) = -^ J dt' J dt" r Ti exc [V(t'),[V(t"), p(t)]], 
t t 

(A3) 

where the collision time r c is of the order of the scatter- 
ing length divided by the thermal velocity and is some 
nanoseconds in duration. The excited states, {nk}, are 
not much influenced by a single collision and the density 
matrix p(t) is taken to be p(t) = p YQS (t) (g> p(&o,aJ,i). It 
is important to note that the trap frequencies are of or- 
der a few hundred Hertz and the integrals over t' and t" 
can be simply replaced by The rate of collisions r 

is governed by the particle mean free path, the thermal 
velocity and the number of (density of) excited states. 
The rate r will contribute only to the constant k of e.g. 
Eq. (|A4|) and is not of interest here. 

Using Eq. l)A3jl and the interaction Hamiltonian Eq. 
(|14fl we obtain the following coarse grained equation of 
motion for the reduced density operator of the interacting 
Bose gas 

p(a , a , t) = rSp = --^] Tr{/3k/3^p,. es }[a k a k aodJpW- 
a k alp(t)a al + adj.]- 



- ^ Tr{/3 k /3 k pr es }[aka k 4ao/o(£)-ak«op(t)4ak + adj.]. 

k^O 

The constant k is an overall rate which plays no role in 
the present problem. 

In this paper, we are most interested in the atomic 
statistics of the condensate and hence define the joint 
probability for having no atoms in the condensate and 
m,ri2,.., rtk,-- = { n k} atoms in the various excited states 
as: 

P({n k },n ) = ({n k },no|p|{n k },no) (A5) 
In view of the definition l|A5(l . Eq. 1A4II implies 

P({n k };no) = -|^((/3 k /3k) + l)[n k K+l)P({n k };no)- 

k^O 
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nkn -P({n k }; n Q - 1) + adj.]- 



Defining 



\ £ <ALAc>[K + l)n P({n k }; n ) 

k#0 



(n k + l)(n + l)P({n k }; n + 1) + adj.] (A6) 

The condensate statistical distribution P(no) is then 
related to the joint probability distribution P({n k };n ) 
and the conditional probability P({n^}\nQ) vi a Bayes 
rule 

P(W}; no) = P({n w }\n )P(n Q ). (A7) 

We obtain P{tiq) from P({n k };no) directly by sum- 
ming over all excited states {n k } 

P(no)= £P({n k };n ), (A8) 
{«k} 

so that we have 

Pn - + l)[(nk)no(«0 + l)Pn - 

k#0 

(n k )„ _in Pno-i]- 



K H^k/3k)[((«k>« +l)noPno-((«k)„o + l + l)K + l)Pno + l] 
k^O 

(A9) 

where we have introduced the convenient conditional av- 
erage notation 



i 



k k/ e phu k _ i ' 
where /3 = l/fc s T. Eqs. (|BTJ) and jE2l 

now read 

k/0 



and 



H no = £?7fc(l + (n k ) no ). 

k^O 



(B3) 



(B4) 



(B5) 



We then approximate the conditional thermal average as 



(nk)n Q = (N - n )- 



n k 



where the usual atomic thermal average is given by 

1 



nk 



(B6) 



(B7) 



Now an ensemble of atoms and phonons in an ordinary 
gas obeys the simple rate equation 



dt 



n k = ~f(n k + l)i] k - jn k (r) k + I), (B8) 



where 7 is an uninteresting rate factor. At steady state 
Eq. (EHJ yields 



(n k + l)r] k = UkiVk + 1), 



(B9) 



£ra k P({ra k }|ra ). 



(A10) 



which implies that r\ k = n k and we may write Eq. 1|B6|) 

as 



(n-fc)no = (N - no) 



Vk 



The master equation obtained here is diagonal. More 
general master equation can be off-diagonal, as we discuss 
in Appendix C. 



APPENDIX B: THE CROSS COUPLING 
COEFFICIENTS 

In this appendix, we present a brief reminder of how 
the heating and cooling coefficients H no and K no were 
handled in CNB2 for an ideal gas. There, we encountered 
the heating and cooling coefficients 

#no = £K> no (l + (S k S k )), (Bl) 



J2k^O Vk 



(BIO) 



Hence we have 



Y t < n ^Vk = ^ 1 ^ y Z4 = (N-no)v (Bll) 



TV 

k^O k^O 

where we have introduced the notations 



H = £ f] k and 



V 



k^O 



H ' 



(B12) 



In this way we arrive at the K no and H no coefficients of 
CBN2 given by 



and 



k^o 



#»o=E(^k)(l + (n k )„J 

k#0 



and 



(B2) 



K no = (N-n )(l + r)), 



H no =H + (N- no)v- 



(B13) 



(B14) 
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APPENDIX C: MASTER EQUATION 
GENERALIZATIONS 



-[Hi + (N - n )r]i}nop no . 



Our interaction Hamiltonian can yield an off-diagonal 
master equation. That is, the main working master equa- 
tion in operator form is 



1 dp(t) 
k dt 



+ [Hi + {N -n -l) m ](n + l)p 



n + l,rao + l" 



= - < 4ak/3k/3£ >n {a alp{t)-al P (t)a )+adi. _ j \^/n Q (n - l)p„ - 2 ,n + + 2)(n + lK ,„ + 2 



X! < PtPkO, k al > na (aja /5(t) - a /3(i)aj) + adj. 



k#0 



-2\/no(no + l)p„ _i, no+ i 



^2 < /3k/3l k "ka-k >n {%alp(t) - aj/3(t)aj) + adj. 



k#0 



\/ (no + l)(no + 2)p„ 0+2 ,no + V n o( n a - l)p n ,n — 2 



- ^ < Ac/3-ka k ai k (aoaop(t) - a p(t)a ) + adj. 

k#0 

(CI) 

New contribution of the form ajajpo^) — ajp(i)aj and 
doaQp(t) — dop(t)ao yields off-diagonal terms in the den- 
sity matrix equation. Taking into account that 



< $kA-k >=< «ka-k >= u k v k 



we obtain a master equation which couples the diagonal 
and off-diagonal terms in the density matrix 



1 

2k' 



— Pno-na = ~(N - Uq)(1 + T]l)(n + l)pn ,n + 



+ {N - n Q + !)(! + Vi) n oPno-i,n a -i- 



-2i/no(no + l)p na +i 



no — 1 



(C2) 



where TLi and rji are given by previous formulas (|29|l and 
(JSTJ, and 



J 



E2 2 

k/O 



1 



1 



(C3) 



The more general analysis of fluctuations with off- 
diagonal terms will be presented elsewhere. We have also 
studied a two atom master equation in which two atoms 
at a time are added or removed from the condensate. 
This is analogous to the two photon laser. The results of 
such an analysis are in basic agreement with the present 
findings and will be published elsewhere. 
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